Apparent electron-phonon interaction in strongly correlated systems 
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We study the interaction of electrons with phonons in strongly correlated solids, having high-Tc 
cuprates in mind. Using sum-rules, we show that the apparent strength of this interaction strongly 
depends on the property studied. If the solid has a small fraction (doping) 5 of charge carriers, 
the influence of the interaction on the phonon self-energy is reduced by a factor 5, while there 
is no corresponding reduction of the coupling seen in the electron self-energy. This supports the 
interpretation of recent photoemission experiments, assuming a strong coupling to phonons. 
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There has been much interest in electron-phonon 
coupling in high- Tc cuprates after photoemission spec- 
troscopy (PES) studies of Lanzara et al. showed a 
strong coupling to a mode at 70 meV. This was inter- 
preted as a coupling to a half-breathing phonon, where 
O atoms in the Cu02 plane move towards Cu atoms. 
This interpretation is supported by the anomalous soft- 
ening under doping ^, ^ A 111 and large width of this 
phonon. It is intriguing, however, that a substantially 
larger apparent coupling strength A was estimated from a 
kink in the PES spectrum than what is suggested from 
the phonon width and softening 0. Similar considera- 
tions apply to the Big phonon width and the anomaly 
at 40 meV seen in PES 0. These estimates of A were 
based on theories which assume noninteracting electrons. 
Using sum-rules, we show that for a strongly correlated 
system the apparent A deduced from such theories should 
depend strongly on the property studied. We consider a 
doped Mott (charge transfer) insulator, such as a cuprate, 
which has a fraction S carriers (where S typically is small) . 
The influence of the electron-phonon interaction on the 
phonon self-energy, determining its width and softening, 
is then reduced by a factor of the order of 5 compared to 
a system without electron-electron interaction. For the 
electron self-energy, determining the photoemission spec- 
trum, there is no comparable reduction. This explains 
why the A deduced from the phonon self-energy appears 
smaller than the one deduced from PES, and it supports 
the scenario that phonons give a large contribution to 
structures in PES. For similar reasons, there should be 
no reduction ~ (5 in the phonon induced carrier-carrier 
interaction. To address these issues, we present a method 
for determining A in exact diagonalization approaches. 

The electron-electron interaction can strongly reduce 
the effects of the electron-phonon interaction on phonons 
. The electron density is rearranged in response to the 
excitation of a phonon, and this rearrangement acts back 
on the phonon, contributing to the width and energy shift 
of the phonon. The response of the electrons depends on 
electron hopping, which is hindered by interaction effects. 
In particular, if the interaction is strong enough to lead 
to a Mott insulator (for (5 = 0), the electron-phonon con- 
tribution to the phonon width vanishes. 

For cuprates, this can be studied in the t-J model [l(Tj |. 



which has one site per Cu atom. Each site is occupied by 
either a Cu Sd-hole or a Zhang-Rice singlet, composed 
of a Cu 3d-hole and an O 2p-hole. The doping 6 gives 
the fraction of singlets, which provide the charge carriers. 
The Haniiltonian is given by 
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where c]^ creates a d-hole on site i if this site previously 
had no hole. A hole can hop with the hopping integral 
t to sites occupied by singlets and vice versa. The spins 
of the 3d-holes have a Heisenberg interaction with the 
interaction strength J. In the t-J model with phonons 
[ll. 12. 13, .14, 15, 16} , the phonons couple mainly to the 
on-site energies and only weakly to the terms describing 
hopping and spin-spin interaction |llL Il5| | . In the follow- 
ing we only include the coupling to the on-site term, 

= E 5q("^ - l)(&q + bl^y'' '''^ (2) 



where N is the number of sites, is a coupling constant. 
Hi measures the d-hole occupation on the site at and 
6q annihilates a phonon with a wave vector q. 

The phonons only couple to sites with Zhang-Rice sin- 
glets, i.e., sites with no d-holes. The phonon self-energy 
n(q, (jj) can then be expressed in terms of the charge- 
charge response function xili^)- 
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where Do{c{,i-u) is the free phonon Green's function. 
KhaliuUin and Horsch ^3] showed that there is a sum- 
rule 
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|Imx(q,w)|dw = 25(1 - (5)7V, (4) 



This result is a factor 2(5(1 — S) times the result for non- 
interacting electrons in a half-filled band. Since x(q, a;) 
becomes small for small (5, the same is true for n(q, a;). 
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In this hmit the denominator in Eq. (PJ is not very im- 
portant, and the sum-rule in Eq. Q also applies approx- 
imately to n(q, a;)/(7q. 

To understand the result in Eq. Q , we notice that the 
system can respond to the perturbation of a phonon by 
transferring singlets to sites with Cu Sd-holes. If there are 
few singlets, i.e., S is small, the response of the system is 
weak and the phonon self-energy is small. Since typically 
S ~ 0.1, this drastically reduces the phonon softening and 
width. 

It is interesting to study PES and inverse PES (IPES) 
to see if many-body effects also drastically change the ef- 
fects of the electron-phonon interaction in these cases. 
These spectra are described by the spectral function 
A(k, w) = ImG(k, w — i0"'")/7r, where k is a momentum 
and uj a frequency. G(k, z) is the electron Green's func- 
tion 



G(k, z) = 
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where ak is a weight, z is a complex number and S(k, z) 
is the electron self-energy. The z-independent part of E 
is included in the energy Ski so that S(k, z) ~ b]s^/ z for 
large z. To determine 6k, we expand the Green's function 
in 1/z 
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where (w")k = / w"A(k,w)du;/ / yl(k,w)dw. The coef- 
ficient &k can then be expressed in terms of (a')k and 
(a-'^)k- Using the analytical properties, we can relate 6k 
to a sum-rule over Im I](k, a; — iO^). Since there are no 
phonons in the ground-state in the model Q for (5 = 0, 
we can easily calculate expectation values of powers of 
ifep, needed to obtain (a;")k- From this we deduce a 
sum-rule for the difference Egp between the self-energies 
with and without electron-phonon interaction. 
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which is valid for (5 = 0. The integration only runs over 
energies corresponding to the PES spectrum, since the 
weight of the IPES spectrum is zero for 5 — 0. The 
result lO is equal to the lowest order result (in g^) for 
nonintcracting electrons. While the effect of the electron- 
phonon coupling is strongly reduced by the small doping 
for the phonon self-energy, there is no such reduction for 
the electron self-energy. Re Eep is related to Im Sgj, via 
a Hilbert transform. 

The sum-rule in Eq. (UJ can be understood if we no- 
tice that a singlet with the wave vector k, created in PES, 
can easily be scattered by phonons to other states k -I- q. 



since only a fraction (1 — 5) of these are occupied by other 
singlets. We therefore expect a strong electron-phonon 
interaction in PES also for a finite but small S. As a 
result, one would then expect that the electron-phonon 
interaction appears to be a factor of l/(c(S) stronger in 
Re E(q, w) than for the phonon width 2 Im n(q, w), 
where c ~ 2 — 4 depends on the assumptions about the 
tj-dependencies of Im E and Im 11. Although the argu- 
ments above show that the right hand side of the sum- 
rule in Eq. Q should not go to for (5 ^ 0, the result 
is, nevertheless, highly nontrivial. The right hand side is 
independent of k, t and J. It is also interesting that it 
remains proportional to for large g. 

From these arguments it follows that the phonon in- 
duced attractive carrier-carrier interaction should also 
be effective, since for small 6 the carriers (singlets) can 
scatter each other via phonons with few restrictions. 
This may be helpful for superconductivity. In partic- 
ular, it should be possible to have a strong phonon in- 
duced carrier-carrier interaction, without the correspond- 
ing phonon going soft. 

To illustrate these points, we study the t-J model with 
phonons, using exact diagonalization. We consider a fi- 
nite cluster with 4x3 sites, and include the entire breath- 
ing phonon branch. The Hilbert space is limited by only 
allowing for states which have a maximum of K phonons, 
where typically K — 7. Ideally, we would calculate the 
PES spectrum as a function of k. We would then ex- 
pect to see a kink in the dispersion at energies of the 
order of ujph away from the Fermi energy, where ujph is a 
typical phonon energy. For the small clusters which can 
be studied in exact diagonalization, this approach is not 
possible, due to the few k-points available. 

Instead we focus directly on the self-energy. The spec- 
tral function A(k, oj) is calculated using exact diagonal- 
ization. The Green's function is then obtained from a 
Kramers-Kronig transformation. Finally, E(k, w) is ob- 
tained by inverting Eq. (jSJ. In a similar way we deter- 
mine n(q, w) from the phonon spectral function. The w- 
dependence of Re E(q, w) then gives information about 
the kink in PES and Im n(q, cjp^) gives the phonon 
width. 

Due to the small size of the cluster there are few many- 
electron states in the phonon energy range, which makes 
it hard to extract phonon widths, even in the approach 
above. The main points of the paper can, however, also 
be illustrated by using a larger phonon energy. We have 
therefore increased the bare phonon energy to uj^^f^ = 0.5 
eV. This requires a corresponding increase in the cou- 
pling constant. We have chosen a multiplying factor 
of 3.4, which leads to an apparent coupling st reng th to 
the phonons of the order seen experimentally We 
have furthermore used the parameters t — OAT eV and 
J/t = 0.3. 

Figure^shows the frequency integrals of the imaginary 
parts of the electron and phonon self-energies as a func- 
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FIG. 1: Frequency integrals over the imaginary parts of the 
phonon self-energy for q — 7r/(2a)(l,0) and the electron- 
phonon contribution to the k-averaged electron self-energy. 
The results are divided by the coupling constants and they 
were obtained for a 4 x 3 cluster with periodic boundary con- 
ditions, n was calculated for S — 1/12 and Sep for 5 = 0. The 
self-energies were given a 0.2 eV (FWHM) Lorentzian broad- 
ening. The figure illustrates how these quantities converge to 
approximately 25 ^ 0.17 [Eq. @] and unity [Eq. 0] for the 
phonon and electron self-energy, respectively. 



tion of the upper limit lo. The limit — )■ oo corresponds 
to the sum-rules (0)) and {Tj). To obtain dimensionless 
quantities, the coupling constants have been divided out. 
The figure illustrates the large ratio, ~ 1/(2(5), of the 
sum-rules. The sum-rule in Eq. applies to an aver- 
age over q 7^ of X, but we found that it is also rather 
accurate for U/g"^ for an individual value of q, as is illus- 
trated in Fig. ^ The reason is that the denominator in 
Eq. ||2J) is not very important for the coupling strengths 
used here and that the sum-rule for an individual q is 
not very different from the average over q. 

The coupling strength A is often determined from the 
phonon width 7 = 2 Im Il{q,ujph) 0| or its softening, 
Aujph = Re Il{q,ujph), using 
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appropriate for noninteracting electrons. Here N{0) is 
the electron density of states per spin and a 1 depends 
on the precise w dependence of Im n(q, w). Figure [21 
shows the broadened Im n(q, uj). The broadening (0.4 eV 
FWHM) was chosen in such a way that the fine structures 
due to the finite cluster size were removed and so that 
the expected behavior Im H(q, cj) ~ to for small lu was 
obtained. The phonon is softened to ujph — OA eV, due 
to the electron-phonon interaction. For this frequency we 
obtain the FWHM of the phonon as 7 = 0.08 eV. This re- 
sult depends on the broadening of Im H(q, oj) and results 
differing by ±30% could be obtained for other reasonable 
broadenings. Based on the width of Im H(q, w), we esti- 
mate N{0) ~ 0.5 states per eV and spin and A'^ = 0.2. 




FIG. 2: Im n(q,a;) of a 4 x 3 cluster for q = 7r/(2a)(l,0) 
and S = 1/12. The self-energy has been given a Lorentzian 
broadening of 0.4 eV. 



From the phonon softening Aujph/ujpii — 0.2, we obtain 
a similar result for A'^. 

We emphasize the difference in this approach, using a 
broadened H, from an approach where the phonon spec- 
tral function is broadened until a smooth spectrum is ob- 
tained. The latter approach leads to an additional width 
of the peaks. For the present system, such a large broad- 
ening would be required that the approach would not be 
meaningful. The broadening of Im H(q, w), on the other 
hand, does not generally add to the width of the phonon 
spectral function, since it essentially only distributes the 
contributions to Im H(q, a;) more uniformly in co. 

We next consider the coupling strength seen in the 
electron self-energy determining PES. We can determine 
a A^ as 



A^ 



dReI]ep(w) 



dui 



j=0 
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This leads to A^ = 0.6, which is about a factor of l/{cS) 
larger than A'^ with c = 4. As can be seen from Fig. 
Im E is small for |w| < ojph = 0.4 eV. For this frequency 
range we find that Re S(k, w) ~ —X^ut. 

To summarize, we have found that in strongly corre- 
lated systems the apparent strength A of the electron- 
phonon coupling crucially depends on the property of 
interest. For the t-J model of a high-Tc cuprate with a 
fraction S carriers, sum-rules for the imaginary parts of 
the electron and phonon self-energies show a reduction 
by a factor S for the phonon but not the electron case. 
This suggests that the apparent A deduced from phonon 
widths and softenings is reduced by such a factor, while 
there is no reduction in the electron self-energy. This pro- 
vides support for phonons being essential for kinks seen 
in photoemission. Similar arguments suggest that the 
phonon induced interaction between the carriers is not 
reduced by a factor of 6, which may be of importance for 
superconductivity. 
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